On the indefinite Helmholtz equation: Complex stretched absorbing boundary layers, iterative analysis, and preconditioning
نویسندگان
چکیده
This paper studies and analyzes a preconditioned Krylov solver for Helmholtz problems that are formulated with absorbing boundary layers based on complex coordinate stretching. The preconditioner problem is a Helmholtz problem where not only the coordinates in the absorbing layer have an imaginary part, but also the coordinates in the interior region. This results into a preconditioner problem that is invertible with a multigrid cycle. We give a numerical analysis based on the eigenvalues and evaluate the performance with several numerical experiments. The method is an alternative to the complex shifted Laplacian and it gives a comparable performance for the studied model problems.
منابع مشابه
A new level-dependent coarse grid correction scheme for indefinite Helmholtz problems
Pushed by the rising interest in high resolution requirements and high-dimensional applications, the diffusion term in the Laplacian equation drives the condition number of the associated discretized operator to undesirable sizes for standard iterative methods to converge rapidly. In addition, for realistic values of the wavenumber k(x) in (1), the Helmholtz operator H becomes indefinite, destr...
متن کاملA preconditioned iterative solver for the scattering solutions of the Schrödinger equation
The Schrödinger equation defines the dynamics of quantum particles which has been an area of unabated interest in physics. We demonstrate how simple transformations of the Schrödinger equation leads to a coupled linear system, whereby each diagonal block is a high frequency Helmholtz problem. Based on this model, we derive indefinite Helmholtz model problems with strongly varying wavenumbers. W...
متن کاملGMRES-based multigrid for the complex scaled preconditoner for the indefinite Helmholtz equation
Multigrid preconditioners and solvers for the indefinite Helmholtz equation suffer from non-stability of the stationary smoothers due to the indefinite spectrum of the operator. In this paper we explore GMRES as a replacement for the stationary smoothers of the standard multigrid method. This results in a robust and efficient solver for a complex shifted or stretched Helmholtz problem that can ...
متن کاملCompressed Absorbing Boundary Conditions via Matrix Probing
Absorbing layers are sometimes required to be impractically thick in order to offer an accurate approximation of an absorbing boundary condition for the Helmholtz equation in a heterogeneous medium. It is always possible to reduce an absorbing layer to an operator at the boundary by layer-stripping elimination of the exterior unknowns, but the linear algebra involved is costly. We propose to by...
متن کاملPreconditioning Helmholtz linear systems
Linear systems which originate from the simulation of wave propagation phenomena can be very difficult to solve by iterative methods. These systems are typically complex valued and they tend to be highly indefinite, which renders the standard ILU-based preconditioners ineffective. This paper presents a study of ways to enhance standard preconditioners by altering the diagonal by imaginary shift...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- J. Comput. Physics
دوره 229 شماره
صفحات -
تاریخ انتشار 2010